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SUMMARY 

By the use of orthogonal polynomials developed for discrete sets of data, 
the least-squares equations for determining the optimized stress-rupture para- 
metric constants are obtained in nearly uncoupled form; thus the use of high- 
degree polynominals is pennitted without the loss of significant figures. 
Optimum values of the constants can thereby be accurately obtained. The method 
is applied to the data obtained from the German cooperative long-time creep 
program by using a general parameter of which the Manson-Haferd and Larson- 
Miller parameters are special cases. Good correlation was obtained. An analy- 
sis is also made of creep data obtained for columbium alloy FS-85 with good 
results. A complete Fortran IV computer program is included to aid those wish- 
ing to use the method. 


INTRODUCTION 

One method of extrapolating short-time creep-rupture data to predict long- 
time life involves the use of a time -temperature parameter. This concept is 
based on the assumption that all creep-rupture data for a given material can be 
correlated to produce a single "master curve" wherein the stress (or log stress) 
is plotted against a parameter involving a combination of time and temperature. 
Extrapolation to long times can then be obtained from this master curve, which 
can presumably be constructed by using only short-time data. Three well-known 
parametric methods are the Larson-Miller, Manson-Haferd, and Dorn parameters 
(refs. 1 to 3). These parametric methods have the great advantage, at least in 
theory, of requiring only a relatively small amount of data to establish the 
required master curve. 

More recently a general creep- rupture parameter was introduced by one of 
the authors (ref. 4) that includes most of the currently used parameters as 
special cases. The analysis in the present paper is therefore based on this 
general parameter. 



A significant advance in the practical application of the parametric 
methods was the development of an objective least-squares method for determin- 
ing the optimum values of the parametric constants without plotting and cross - 
plotting the data and without the use of judgment on the part of the analyst 
(ref. 5). This least-squares method involves^ however^ several practical dif- 
ficulties that arise from the fact that in fitting the master curve by a poly- 
nomial^ the set of linear algebraic equations for the coefficients (the normal 
equations) are very ill-conditioned. The determinant of these equations can be 
shown to be related to the Hilbert determinant (ref. 6)^ which rapidly ap- 
proaches zero as its order increases. Thus for polynomials above the second 
degree^ it is necessary to use double-precision arithmetic (16 significant 
digits or more) on the computer^ and for the fifth degree and above the results 
become uncertain even with double-precision arithmetic. This difficulty is 
inherent in the normal least- squares equations and is not limited only to the 
stress -rupture problem. 

The present report presents a method for avoiding the above difficulty by 
using orthogonal polynomials in the representation of the master curve (appen- 
dix a). The use of orthogonal polynomials for representing discrete sets of 
unequally spaced data is described in reference 6 and in more detail in refer- 
ence 7. A further improvement can be obtained by perfozming a linear trans- 
formation on the stresses (or the logs of the stresses) so that all the values 
of stress (or log stress) lie between 2 and -2^ as recommended in reference 7. 

As a result of these innovations^ it became possible to perform all the compu- 
tations in single-precision arithmetic (eight significant digits) up to 18th 
degree polynomials without appreciable round-off error. 

In addition^ this report contains a complete analysis^ in which the gen- 
eral parameter was used^ of all the data for three steels that were obtained by 
NASA through the cooperation of Dr. K. Richard of Faberwerke Hoechst in 
Frankfurt and that were investigated in a long-time cooperative creep program 
in Germany. Some of the data from the latter investigation are included in this 
paper. 

Finally it is shown by means of a concrete example how the parameter tech- 
niques can be applied to creep data to predict long-time creep. For this pur- 
pose the data for columbium alloy FS-85^ as reported in reference 8^ are used. 

A complete Fortran IV program^ as used on the IBM 7094 computer in making 
the calculations^ is presented in appendix B. This program can be used for the 
objective analysis of any set of creep-rupture data by the Larson-Miller^ 
Manson-Haferdj, or the more general parameter of reference 4. 


SYMBOLS 

A^B linear transformation coefficients 

a_,b,c elements of coefficient matrix 
D s t and ar d deviation 
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K degree of freedom 

m degree of polynomial 

n mmiber of data points 

P(a) creep-rupture parameter 
0, polynomial 

q stress exponent 

r temperature exponent 

S sum of squares of residuals 

T temperature 

temperature intercept 
t time to rupture 

ta time intercept 

u coefficient of polynomial function 

X scaled log stress 

X log stress 

y log time 

log time intercept 
constants from recurrence relation 
a stress 

T a<l(T - Ta)'' 

Subscripts: 
max maximum 

min minimum 


PROCEDURE 
Ceneral Parameter 


The general creep-rupture parameter introduced in reference 4 has the fol 


lowing form 


i2S-l . log 

gq. 

(T - 


( 1 ) 


where log q, and r are material constants to he determined from the 

available experimental data. The parameter P(cr) is a function of the stress 
and, when plotted against stress, is referred to as a master curve (fig. 1^ 
p. 9). If q = 0 and r = 1, the Manson-Haferd parameter is obtained. If 

q = 0, r = -1, and = -460^ F, the Larson-Miller parameter results. If 

q = 1 and r = 1, the stress-modified parameter suggested in reference 9 is ob- 

tained. Finally, if q = 0, equation (l) reduces to the parameter proposed by 
Manson and Brown (ref. 10). 


The object is to find the best values of the constants q, log 
and r so that the master curve best fits the data. To find these values, the 
method of least squares is used whereby the master curve is represented by a 
polynomial in the logarithm of the stress, and the best fit is obtained by mini- 
mizing the sum of the squares of the deviations (the residuals) of the data 
from the curve. The calculation procedure will now be described. The details 
of the derivation are given in appendix A, and a Fortran IV computer program 
using this method is given in appendix B. 


Calculation Procedure 


To simplify the notation, the following symbols are introduced: 


T = a<l(T - Ta)^ 
y = log t 
X = log a 

7a = “ta 


Then from equation (l) it follows that 

y = crSa + tQ(x) 


> ( 2 ) 




(3) 


where in reference 5^ Q(x) was represented by a simple polynomial of the form 


Q(x) = sq + aj_x + agx^ + . . . + a^jX™ (4) 

The least-squares equations obtained sometimes led to difficulties as indicated 
in the IITORODUCTIOH. These difficulties can be avoided, however, by rewriting 
equation (4) in terms of polynomials that are orthogonal over the set of data, 
as defined in appendix A. Thus assume 
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m+1 

Q(x) = u^Q^(x) + UgQgCx) + . . . + ^ u.Q.(x) (5) 

0=1 

where is an unloiown constant, m is the degree of the highest degree 

polynomial, and Qj(x) is a polynomial of degree j - 1 that satisfies the 
orthogonality conditions described in appendix A. The use of orthogonal poly- 
nomials permits the solution of the least-squares equations directly in closed 
form, thus the loss of a large number of significant digits is avoided. The 
method of calcixLating will be discussed in appendix A. 

If equation (5) is substituted into equation (s)^ an equation with m + 5 
unknown constants results for the case of the general parameter. For the case 
of the linear parameter there are m + 3 constants^ and for the Larson-Miller 
parameter there are m + 2. It is necessary that the number of data points n 
always equals or exceeds the nutliber of unknown constants. 

The constants are determined so that equation (3) fits the data best in 
the least-squares sense. To accomplish this^ the sum of the squares of the 
deviations is minimized; that is^ 

® ^ E [^i ■ ‘"K - ^i^(^i)] (6) 

i=l 

is made a minimum. Because the equations are nonlinear in some of the unknown 
constants a trial and error procedure must be used. A set of values is assumed 
for q^ r^ and Ta. and the corresponding best values of yg, and uj are de- 
termined. A different set of values for q, r, and Tg is then chosen, and 
again the best values of y^ and uj are calculated. Several sets of values 
of q, r, and Tg are tried, and the values corresponding to the overall best 
fit are determined. For the case of the linear parameter, only the value of 
Tg, is varied (q is always equal to zero, and r is always equal to l). For 
the Larson-Miller parameter, Tg is equal to -460° F, and no trial and error 
procedure is needed. 

As a measure of the fit, the standard deviation D, defined by 
is used^ where K equals 


m + 5 

general parameter 




m t 3 

linear parameter 


► 

(8) 

m + 2 

Larson-Miller parameter 

> 




The smallest value of D will correspond to the best fit. 
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To determine the best values of and uj for a given set of values of 

‘ttie following calculations are made- Firsts the logarithms of 
the stresses are scaled so that they lie in the range -2 to 2^ as suggested in 
reference 7. The reason for this is discussed in appendix A- Thus define a 
variable X by 


A = 


X = Ax + B 
4 

^ax “ ^in 




B = -2 


^max ^min 
^ax ” ^in 


(9a) 


(9b) 




The polynomials Qj(^) Q'^re now calculated for each of the data points by 
using the following formulas : 


Vi = 

^ X^T?Q^(Xi) 


m > j > 1 


( 10 ) 




i=l 

a. = 

J n 


m > j > 1 




i=l 


n 


V (10a) 




I'} - 


i=l 


i=l 


m>j>l^ Q-, =1^ and (3-, = 0 


where n is the number of data points^ is the scaled value of log for the 
i*^^ data pointy and Tj_ is equal to ^^(^q “ data point for 

the chosen values of Tg^^ q^ and r. 

It is to be noted that the degree of the polynomial Q(x) of equation (5) 
can be increased by merely computing the next polynomial in the series Qj^4.2 

without having to recompute any of the previous ones. This is one of the ad- 
vantages of using orthogonal polynomials. 

Once the values of Qj have been computed for each of the data points^ 

y„ and u^ can be calculated as follows: 
a j 
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I 



Let 


n 



i=l 


n 


a . 

3 


= 2 


i=l 


n 


- E 


i=l 


=0 = X 

i=l 


vhere j = 1^2 
Then 


m + !• 


n 




i=:l 


m+1 


y« = 


, _ y 

'0 Aj b. 
j=i J 


m+l p 
J=1 


J ~ b. 




> (11) 




y ( 12 ) 


Note that if q = 0^ equals the number of data points n. Thus by means of 

equations (9) to (l2)^ the best values of ya and Uj to fit the data are 
found for a given choice of T^^ q^ and r. The Fortran IV program described 
in appendix B automatically scans all the desired values of Ta^ q.^ and r and 
chooses the best set from all the submitted values as determined by the smal- 
lest value of the standard deviation D, as defined by equation (7). The 
method can be illustrated by a simple example: consider a set of theoretical 

data^ vhich fit the following equation exactly 
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( 13 ) 


9^5 - log t ^ 10-3(7.02 + 0.467 x + 0.061 x^ + 0.00928 x3) 

T - 600 

Eight data points satisfying this equation are given in columns 2 to 6 of 
table I. For this data ^ 1°S "^a = ya = 9*3> Suppose, however, 

that these eight data points were obtained experimentally and that the values 
of Ta log tg^ were not known. The problem then is to find the best 

values of Ta and log t^ to fit the data by the linear parameter. These 
values can readily be found by using the equations of the previous section. 
First, from column 6 of table I 

(log a)j^g^ = 4. 75051 
(log = 1.81954 


Therefore from equations (9b) 


A = 1. 36474 
B = -4.48319 

and by means of equation (9a) the were computed and are given in column 8. 

For illustrative purposes three values of Ta were chosen, 500°, 600°, 
and 700° F. For each of these values of Ta, values of Ti, aj, Pj, and 
Qj(% ) were computed by means of equations (2), (lO), and (lOa), and the values 
of aj, bj, and cj were computed by equations (ll). The results are tabu- 
lated for Ta = 600° in columns 9 to 12 of table I and in table II up to a 
third degree polynomial. 

The values of ya and uj were then computed by using equations (12) for 
each of these three values of Tg^ by first assuming m = 2, then m = 3, and 
finally m = 4, corresponding to polynomials of second, third, and fourth de- 
grees, respectively. For each of these cases the standard deviation D was 
computed from equation (7) with S being given by equation (6) and Q by 
equation (5). The results are summarized in table III. The least value of D, 
signifying the best fit, is obtained for m = 3 and Ta = 600° F. The cor- 
responding value of ■ yg^ is 9.5. These values, of course, correspond to equa- 
tion (13), from which the data were generated. 


Application to Data from German Cooperative Long-Time Creep Program 

As part of the German cooperative long-time creep program, a sufficient 
amount of material of each of three steels was supplied to NASA to permit the 
running of short-time tests necessary to predict the results at long times ob- 
tained in the German test program. The composition of these steels is shown in 
table IV. 

The resvilts of the NASA tests, which were used in the subsequent analysis. 
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Figure 1. - Master curve for steel K (27b KK), calculated 
from NASA data between 10 and 3700 hours. 


are slio-wn in table V. Table VI shows the 
results of the long-time German test pro- 
gram. The three steels will be designated 
briefly as steel K, steel and steel P. 

With the use of the test data shown 
in table V a complete analysis was made 
by the previously described method. The 
general parameter discussed in the 
INTRODUCTION was used^ and the best values 
were obtained for the parametric constants 
for each of the three steels. 

All the data obtained for these 
steels are shown in tables V and VI. 

Many of the data points were obtained for 
purposes other than the application to 
time-temperature parameters^ as described 
in this report. As already discussed in 
references 4 and 11^ a much smaller 
amount of data is needed when an accel- 
erated program is desired; however^ since 
these data were already available^ all the 
data indicated in tables V and VI were 
used to obtain the best possible para- 
metric constants. 


For all three steels the analysis showed the stress exponent 1 to be 
zero^ but the temperature exponent r to be different for each of the three 
materials. For steel K the best value of r was 1^ which indicated that the 
best fit is obtained by the linear parameter. For steel P a value of r of 
-1 was obtained^ which indicated a parameter similar to the Larson-Miller param- 
eter; however^ the corresponding value of Ta was 200^ F rather than -460^ F 
used in the Larson-Miller parameter. For steel C the value of R was 2.5. 


Figure 1 shows the results for steel K. Here the master curve consists of 
a plot of stress against the optimized parameter (T - 300)/(log t - 16.54). 

Figure 2 shows the isothermals computed by using the optimized parameters^ 
as shown on each of the figures. The range of the NASA data used to obtain 
these parameters is also shown on each of the figures. The data points shown 
are the German results obtained to date. The predictions up to 100 000 hours 
from the NASA data based on the optimized parameters agree well with the 
(German data, if scatter and differences in testing technique between the two 
organizations are considered. 

Figure 3 shows a comparison for each of the three steels between the best 
linear parameter, the best Larson-Miller parameter, and the best general param- 
eter. Although for some of the steels fair agreement can be obtained with one 
or the other of these parameters, it is clear that the general parameter is 
superior when all the materials are considered jointly. If any one of the 
special cases of this parameter is to be chosen for all materials, the linear 
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Stress, ksi 



(a) Steel K (27b KK); parameter, 
P = (logt-16. 54)/(T-300). 


(b) Steel C (23b CK); parameter, 
P = (logt-8.87)/(T + 100r-5 


(c) Steel P (14a PA); parameter, 
P = (T“ 200)(logt + li. 13). 


Figure 2. - Analysis of German steel data by generalized parameter with optimuip constants (where T is temperature, and t is time to rupture). 



Temperature for 
German data 
above 5 hr, 

o 932 

O 1022 

o 1112 

^ 1202 

□ 1292 



(b) Steel C (23b CK). <0 Steel P (14a PA). 

Figure 3. - Analysis of German steel data by several parameters (where T is temperature, and t is time to rupture). 
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Parameter, P = - 

’ log t - 5. 76 



(a) 5-Percent strain. (b) Master curves obtained for 1- 2- and 5-percent strain. 

Figure 4. - Analysis of creep data for columbium alloy FS-85 by linear parameter. 


parameter would appear to be the best choice. 


Application to Creep Data 

Although there is no fundamental reason why the same parameter is capable 
of representing both creep and rupture data^ it has nevertheless been found 
empirically (refs. 1 and 2) that the dual role of the same parameter leads to 
reasonable results. Experimental data for creep are much more limited^ however^ 
than that for rupture^ and such data tend to contain more scatter; hence ^ 
analysis of creep data by the parametric approach has been limited in the past. 

The method of the present report can be applied directly to creep data 
without any change. All that is necessary is to redefine t as the time to 
attain a specified amount of creep rather than as the rupture time. Thus^ it 
is assumed that for a given amount of creep^ say 1 percent^ a plot of log o 
against a parameter^ such as that given by equation (l)^ will produce a single 
master curve. For a different amount of creep^ say 5 percent^ a different 
master curve can be obtained^ but it is assumed that the parametric constants^ 
such as log tg, and Tg^^ remain the same and that they equal the values ob- 
tained from rupture data. 

Calculations of this type were performed for columbium alloy FS-85. The 
creep tests were limited to runs of approximately 1000 hours; the data are 
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given in table VII^ as tal;en from reference 8. Figure 4(a) shows the data for 
5-percent creep strain^ and figure 4(b) shows the master curves obtained for 
1-^ 2-^ and 5-percent strain as well as the parametric constants obtained by 
the method of this report. While scatter in the creep data is high^ the cor- 
relation must be regarded as good. In general^ the points agree well with the 
master curve. 

Although these results are encouraging^ much more work is necessary before 
it can be concluded that the parametric approach is completely valid for creep 
data. If it is eventually concluded that the parametric approach is valid for 
creep data and in particular that the parametric constants are the same for 
both the creep and rupture processes^ it is obvious that a great saving in test 
facilities and test program planning will result. It therefore seems very 
worthwhile in future studies to give more attention to the correlation and 
extrapolation of creep data by the parametric method. 


Lewis Research Center^ 

Rational Aeronautics and Space Administration^ 
Cleveland;, Ohio^, May 3^ 19G5. 
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APPENDIX A 


ORTHOGONAL POLYNOMIALS AND LEAST- SQUARES DEDERMINATION 
OP PARAMETRIC CONSTANTS 

A set of polynomials Qj (x) are said to be orthogonal over an interval 
with respect to the weighting function t(x) if they satisfy the following 
relation 


rx=Xg 

I T^(x)Qj(x)Qj5.(x)dx = 0 j k (Al) 

*'x=X3_ 

Similarly a set of polynomials can be defined to be orthogonal over a set of n 
discrete points by the following relation 

n 

T?Qj(xi)Qj5.(Xi) =0 J ^ k (A2) 

i=l 

It can be shown (ref. 6)^ that all orthogonal polynomials satisfy a three-term 
recurrence relation of the form 

^k+1 = (^ " “ Pk^k-1 ^ ^ 1 (A3) 

Thus by starting with Qj. = ^ = 0 an infinite set of orthogonal poly- 

nomials can be generated by means of equation (A3) if values for oc^ and 
are known. These can be determined from the orthogonality conditions (eqs. 

(A1) or (A2)). From the relation (A2) it follows that 


n 

^x2Q^(xi)Qk+i(xi) = 0 (Aia) 

i=l 

and 

n 

i=l 

When the recurrence relation (A3) is used to eliminate there is obtained 

n 

2 - PA-i] - 0 

i=l 
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^ - Pksc-l]\-l = ® 

i=l 

When the orthogonality condition (A2) is used^ equations (A5a) and (A5b) reduce 
to 


^ T?(Xi - 0C^)Q| - 0 


1=1 


(A6a) 


1=1 

Solving equations (A6) for cx^ and Pj^ gives 

n 

XiTiQk 


= 0 


1=1 

°^k = — 


z 

1=1 


t2q2 

^l^k 


n 




i=l 


•^k ~ n 


1=1 


2^2 


rQ - 


k-l 


(A6b) 


(A7a) 


(A7b) 


Thus a set of orthogonal polynomials can he generated that are orthogonal over 
a finite set of discrete values of the variable x. Note that these values 
need not be equally spaced^ a condition that is obviously necessary for stress- 
rupture data. 


Scaling of Polynomial Argument 


From the recurrence relation (A3) with Q-i = 1^ it follows that the lead- 


Xf, 


Therefore^ depending on the values of x^^ the 


ing term of is 

values of Qk+i(^i) can become very large or very small. This procedure can 
lead to a loss of significant figures in performing the calculations. It is 
shown in reference 7 ^ by comparison with the Chebyshov polynomials^ that if 
X is scaled so that all the values of lie between 2 and -2^ the polynomial 
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values all le of approximately uniform size. To perform this 

scaling, let ^xaax maximum value of log o and ^±ri ’be the minimum 

value of log a; then let 


X = A log a + B 

(A8) 

2 = -^^^ax ^ 

(A9a) 

-2 = Ax^j^ + B 

(A9b) 


and solving for A and B results in equations (9b). 

It has been found in practice that scaling the values of x as indicated 
does indeed preserve the significance of the calculations. 


Least-Squares Procedure 


In terms of the orthogonal polynomials, equation 


( 3 ) can be written 


m+1 

y = ^ ^ (AlO) 

j=l 

To find the best values of and uj that fit the data, the sum of the 

squares of the residuals is minimized. Thus let 


S = 



a% 



(All) 


Then in order to find the values of and Uj that will mal:e S a minimum, 

S is differentiated in turn with respect to and each u^, and the resulting 

equations are set equal to zero. When this is done, the following set of equa- 
tions is obtained: 


aoya + a^ui + a2U2 + . . 
^l^a + Vl + 0 + . 

0 + ^ 2^2 + • 


+ %+l%l+l = Cq 

. + 0 = 

• + 0 =! Cg 






(A12) 


%+l^a + 0 + 0 + . . . + 


where 
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ao = 


n ^ 



i=l 


n 

aj = ^ aJx^QjCX^) j = 1,2 . . . m + 1 

1=1 

n 

bj = T^Qj(Xj^) j = 1,2 . . • m + 1 

1=1 

°o = X ^^yi 
1=1 

n 

Cj = ^ ^ ^lyi'^j^^l^ J = 1^2 . . . m + 1 

1=1 


>“ (A13) 


J 


It Is to be noted that the only nonzero elements In the coefficient matrix of 
equations (A12) are the diagonal elements and the elements of the first row and 
first column. All the other elements are zero because of the orthogonality 
properties of the polynomials used. This Is one of the major advantages In 
using orthogonal polynomials. In the usual case of data fitting, all the ele- 
ments of the first row and first column, except for the first element, would 
also be zero; and the equations would be completely uncoupled, each uj 
being computed completely Independent of the others, without the necessity of 
solving any sets of equations with the resultant loss of significant figures. 

In this particular case because of the added constant y^, the equations are 
not completely uncoupled, but they are very nearly uncoupled and can readily 
be solved. Thus for any equation after the first 


j ~ b. 


(A14) 


Substituting Into the first equation and solving for ya give Immediately 

m+1 


ya = 


V j j 

^0 - 2 -^ 

j=l ^ 

m+l 2 

j=i J 


(A15) 
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APPENDIX B 


FORTRAN rV PROGRAM 


$ID YAG1202 EKNtST ROlicRTSf JR. - 140 M-S - PAX 6132 
SLIBSIO COMFINUC 

SIBJOb SOURCE 

$IBFTC PRMTRl L 1 S T , KE F , DECK 


c 

LREEP/:>TR6SS~RUPTU^^e PARAMETER PROGRAM 

PRMI 

i 

c 



PKM.T 

2 

c 

NOMENCLATURE IS AS FOLLOWS 

PRMT 

3 

c 



PRMT 


c 

00 

STANDARD DEVIATION 

PRMT 

S 

c 

KK 

DEGREE OF FREEDOM 

PRMT 

0 

c 

KM 

NUMbER OF VALUES OF M READ 

PRMT 

7 

c 

Kvj 

NUMBER OF VALUES OF U READ 

PRMT 

L 

c 

KK 

NUMbER UF VALUES OF R REaU 

PRMT 

) 

c 

KTA 

NUMBER OF VALUES OF TTA READ 

PRMT 

ir 

c 

M 

DEGREE POLYNOMIAL 

PKM1 

1 1 

c 

iM 

NUMBER OF DAIA POINTS 

PRMT 


c 

PP 

PARAMETER 

PRMT 

i 3 

c 

U 

STRESS EXPONENT 

PRMT 


c 

QO 

POLYNOMIAL 

PRMT 

15 

c 

R 

TEMPERATURE EXPONENT 

PRMT 

IG 

c 

RATIO 

ABS( Y-YY)/DD 

PRMT 

i / 

c 

SIGMA 

STRESS 

PRMT 

18 

c 

S IGg 

S 1 oMA* 

PRMT 

iv 

c 

T 

T IME 

PKMT 


c 

TA 

TIME INTERCEPT 

PRMT 

^ 1 

c 

TAU 

SIGMA*»0*{ TT-TTA)**R 

PKMT 

22 

t 

TAUSOR 

TAU*« 2 

PRMT 

^3 

c 

TIME 

calculated T ( iO.*»YY) 

PKMT 

2A 

c 

TT 

TEMPERATURE 

PKMT 

2 S 

c 

TTA 

temperature intercept 

PRMT 

c6 

c 

X 

LOG SIGMA 

PKMT 

^ / 

c 

Y 

LOG T 

PRMT 

2 8 

c 

YA 

LOu TA 

PKMT 

^0 

c 

YY 

CALCULATED LOG T 

PRMT 

3c 

c 



PRMT 

s t 

c 

ALL 

UUANTITIEs in common WITH THIS PROGRAM A\J THIS PAPER 

PRMI 

32 

c 

aTL 

plpresfnted by the same symbol, with repeated 

PRMT 

33 

c 

LETTERS INOILATING IhE UPPER CASE ANU GREEK LETTERS REI... SPLLLL' 

)-PR MT 

3^ 

c 

OUT. 


PRMT 

35 

c 



PKMT 

3 6 

c 

program EXTRAPOLATES CR E E P / S T R E S S-R UP T UR L JATA USING A 

PKMT 

37 

c 

GENERALIZEU PARAMETER 

PRM 1 

38 

c 


PP=(Y/SIGMA^^*0-YA)/(TT-TTA)**R, 

PRMT 

39 

c 

SELECT^ PARAMETER PRODUCING SMALLEST RESIDUAL AND OUTPUTS A 

PRMI 

AG 

c 

COMPLETE 

table, results 0F all other VALUES ARE SUMMARIZED IN 

PKMT 

Ai 

c 

A SHOiTTEK 

table. 

PRMT 

A^i 

t 



PRMT 

A3 

c 

************ I nput ************* 

PRMT 

AA 

c 



PRMT 

A5 

c 

TITLE uAKD, MODE CARD, AND FIVE (5) SETS OF DATA. AT THE ENU OF 

PRMT 

A6 

c 

EACH ScT 

OF data must be a CARU WITH THE WORD *LND* IN THE FIRST 

PRMT 

A7 

c 

IHREE COLUMNS. ALL DATA CARDS (EXCEPTING TITLE AND MODE CARDS) 

PRMT 

A 6 

c 

MUST HAVE 

BLANKS IN THE FIRST THREE COLUMNS. GuLUM,>jS 73-80 ARE 

PRMiT 

A9 

c 

INGORED. 


PRMT 

sG 

c 



PRMI 

:> 1 

c 

TITLE ' ANY ALPHAMERIC INFORMATION — HEADS EACH PAGE OF UUTPuT 

PRMT 

s2 

c 



PRMT 

3 3 

c 

MODE Card 

' ONE OF IHREE WORDS IN COLUMNS i-6, 'LARSON*, 'LINEAR* 

' ,PRMT 

5A 

c 


OR 'GENRAL*. THIS CARD DEFINES 'KK*, THE DEGREE OF 

PRMT 

5 5 

c 


FREEDOM, USED IN CALCULATING GOODNESS OF FIT. 

PRMT 

56 

c 



PRMT 

5 7 

c 

DATA Set 

1 — VALUES UF TTA TO BE INVESTIGATED — ONF PER CARD 

PRMT 

58 
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o o o o 


C 

C 

c 

c 

c 

c 

c 

L 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


L 


FO^Ma'T (3X,F10.0) — 50 VALUES MAXIMUM 

JATA StT 2 — VALUES UF TEMPERATURE EXPO.MENT, R, TJ fiE iNVEbTIGATFU 
UNE PER CAR0**-*FGRMAI { 3 X , F 1 0 . 0 ) '-2C VALUES MAXIMUM 

UATA SeT 3 — VALUES UF STRFSS tXPONItMT,U» TO BE i M V£ T I GA T £ 0 

UNE PER CARD--FORMAT (3X,FL0.0) — 2C VALUES MAXIMUM 

UATA i>tT A-”UEGREES OF POLYNOMIAL, M, TO BE INVESflGATEU 

UNE PER CARO— FORMAT (3X,I2) — MAXIMUM VALUE NOT TO 
EXCEED 20--ZER0 MAY NO f BE USED* 

UATA SuT 5 — UATA POINTS IN THE ORDER TEMPERATURE, STRESS, AND 
IIME--~UisiE :>er PER CARO — FORMAT (3X,3F10*0) 

THE value UF STRESS IS A U TOMA T I C AL L Y OIVIOEU BY lOOC 
FOR ALL CALCULATIONS EXCEPT FINDING THE LOG STRlSS. 
^:00 SET:> MAXIMUM. 

LACH OF THE FIVE SETS UF DATA MUST b6 FULLUWED BY A CARD HAVING 

the word end in the FIRbT ThREE COLUMNS. 

ALL UATA CARuS (EXCEPTING TITLE ANO MOUF CAROS) MUST HAVE THE 

FIRST THREE COLUMNS BLAiNK. 


WITHIN EACH SET, DATA MAY BE IN ANY ORDER. 
IN THE ORDER PRESENTED TO THE MACHINE. 


IT 3E PRUClSSLO 


THE CALCULATIONS ARE PERFORMED IN FOUR (^) LOOPS. 

GOING FROM INNFRMOS) TO OUTFRMOST, THE UJANTITILS ARE VARIED 
IN THE FOLLO/^ING ORDER 

degree POLYNDMIAL, M 
VALUE OF TTA 
TEMPERATURE EXPONENT, R 
STRESS exponent, w 

THE OUTPUT TABLES UTILIZE LFSS THAN 120 COLUMNS ON TriE PRINTER 
AND EXPECT NO CARRIAGE CONTROLS OTHER THAN i, 0, + ANO BLANK. 

A LINE counter IS INCORPORATED TO LIMIT OUTPUT TO 6U LINES PER 
PAGE. FOR EALH N F vV PAGE THE TITLE ANO APPROPRIATE CGLUMM HEADINGS 
ARE PRINTED. PROGRAM EmDS /^ITiI A TRANSFER TO THE INITIAL READ. 

PAGE Counting and error traps most be provided by the operating 

SYSTEM. 

PROGRAM with IBSYS AND lOCsM WILL RUN ON A 16K MACHINE 


PRMT 

59 

PRMT 

60 

PRMT 

61 

PRMT 

62 

PRMT 

63 

PRMT 

64 

PRMT 

65 

PRMT 

66 

PRMT 

67 

PRMT 

6B 

PRMT 

69 

PRMT 

70 

PRMT 

71 

PRMT 

Z2 

PRMT 

73 

PRMT 

74 

PKMT 

75 

PRMT 

76 

PRMT 

77 

PRMT 

76 

PRMT 

79 

PRMT 

I'O 

PRMT 

ni 

PRMT 

o 2 

PRMT 

8 3 

PRMT 

64 

PRMT 

fa5 

PRMT 

66 

PRMT 

87 

PRMT 

bH 

PRMT 

89 

PRMT 

90 

PRMT 

91 

PRMT 

92 

PRMT 

9 3 

PRMT 

94 

PRMT 

9'j 

PRMT 

96 

PRMT 

97 

PKMT 

9b 

PRMT 

99 

PRMT 

lUG 


PRMT 
PRM T 
PRMT 
PRMT 
PRMT 
PRMT 


ICi 
1G2 
1L3 
10^ 
IDS 
i G 6 


LOGICAL TRgGkI, TRGGR2, 1RGGR3 PRMT 

PRMT 

DIMENSION TI1LE(12), tablets, 110), ITBLEI6, 110) PRMT 

PRMT 

EwU I VALENCE ( TA B L E ( i , 1 ) , I T BL E ( 1 , i ) ) PRMT 

PRMT 

COMMON /DATA/SI GMA( 201 ) , T( 201 ) , TT( 20n PRMT 

1 /TRYS/M(21),w(51),R(51),TTA{51) PRMT 

2 /FDATA/SIG'D(200),TAU(200), TAUSURT 200 ) , X { 200 ) , XX ( 200 ) , Y ( 20C ) PRMl 

3 /CALC/PP(2 00),RATI0(200) , TIMt:( 200) , YY( 200) PRMT 

A / t NO/ LND/N/N/DD/ jD/ degree/ DEGREE PRMT 

5 /PLYNML/UT HERK A22 1 ) , Y A , ij T HE R 2 ( 6 3 ) PRMT 

PRMT 

PRMT 

INPUT PRMT 

PRMT 

WRITE (6,999:^) PRMl 

READ (5,9001) ( T I TL L ( K ) , K = I , 1 2 ) PRMT 


I 

lOb 
I 09 

1 A G 
111 
112 
113 
11 ^ 

115 

116 
I 17 
lib 
119 
12»: 
121 
12 ^ 

123 

124 


18 



non oor. no ooo roo r'oo oonoo 


KEAD (5f9001) DEGREE 
K = 0 

10 K = K+1 

READ (!>f9002) CHECK, TTA(K) 

IF (CHECK. i>4E. END) GO TO 10 
KTA = K-i 
K = 0 

15 K = K+1 

READ (5,9002) CHECK, K(K) 

IF (CHECK. NF. END) GO TO 15 
KK = K“1 
K = 0 

20 K = K+i 

READ (5,9002) CHECK, 0(K) 

IF (CHECK. NE. END) GO TO 20 
Kg = K-1 
K = 0 

25 K = K+1 

REAU (5,9003) CHECK, M(K) 

IF ( CHlCK.NF. ENJ) GO TO 25 
KM = K~l 
K = 0 

30 K = K+1 

READ (5,900A) CHf CK , T T ( K ) , S I GMA ( K ) , T ( K ) 

IF ( CHECK. Nt. END) cO TO 30 

N = K-1 

end of input 

FIND LuG stress ANU LOG TIME 

DO 100 K=1,N 

X(K)=AL0Gl0(5lGMA(K) )+3. 

Y( K)=AL0G10( I (K ) ) 

100 CONTINUE 

INITIALIZE constants 

0UI=1. e 5 
LINES=5l 
TKGGR3^.FALSE. 

NTRY=0 

bCALE LOGS OF STRESS 
CALL SCALE 

FIND HIGHEST DEGREE POLYNOMIAL 
MAX = U 

DO 110 K=I,KM 
MAX = MAXO(MAX,M(K) ) 

110 CONTINUE 

MAJOR lUUP - CALCULATES ALL Y(A)S AND RESIDUALS 
WRITES SUMMARY TABLE 
FINDS SMALLEST RESIDUAL 

DO 500 K5=1,KU 

CALCULaTF S IGMA**0 

DO 112 K=1,N 
:,iCG( K) = SlCMA(K) »»w(K5) 

112 CONTINUE 

OU AGO K4=l,KK 


PRMT 

125 

PRMT 

126 

PKMT 

127 

PkMT 

12b 

PRMT 

129 

PKM 1 

13U 

PRMT 

131 

PRMT 

132 

PRMT 

1 33 

PRMT 

1 3A 

PRMT 

135 

PRMT 

1 5t> 

PRMT 

1 37 

P ("v M T 

13b 

PRMT 

13> 

PRMT 

lAC 

PRMT 

lAl 

PRMT 

lA^ 

PRMT 

1 A3 

PRMT 

1 A A 

PRMT 

1A5 

PRMT 

1 A6 

PRMT 

1 AT 

PRMT 

1 A b 

PRMT 

1 A9 

PRMT 

150 

PRMT 

151 

PRMT 

15^ 

PRMT 

15 3 

PRMT 

1 3A 

PKMI 

135 

PRMT 

15o 

PRMT 

l:>7 

PRMT 

I5b 

PRMT 

1 bU 

PRMT 

lt>0 

PRMT 

161 

PRMT 

lo^ 

PRMT 

16 3 

PRMT 

I 6A 

PRM 1 

1 6 •-> 

PRMT 

I6f 

PRMT 

16/ 

PkMT 

1 68 

PRMT 

1 6 9 

PRMT 

1 7C 

PRMT 

171 

PKMT 

1 7^: 

PRMT 

1 73 

PRMT 

1 7A 

PRM.1 

i 7 3 

PRMT 

176 

PRMT 

17/ 

PRMT 

1 78 

PRMT 

1 79 

PRMT 

IbO 

PRMT 

lol 

PRMiT 

18^ 

PRMT 

1C3 

PRMT 

1 8A 

PKMI 

185 

PK^.T 

186 

PRMT 

18 7 

PRMT 

1 h <» 

PRMT 

ib9 

PkN I 

1 -io 
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00 300 K3=1,KTA 

PRMl 191 

c 


PRMT 19^ 

c 

cALCULATe TAU ANo TAU»*2 

PRMT 193 

c 


PRMT l:#4 


UO 120 K=1,N 

PRMT 195 


I JlFF = A6i>{ TT{K)-TTA(K3) ) 

PRMT 196 


IF (TUIFF) 116,115,118 

PRMT 19T 

115 

TAU(K)=0. 

PRMT 198 


OO TO 119 

PRMT 199 

118 

TAU( K) :=SIGQ( K)*T01FF»*,<( K4) 

PRMT 200 

119 

TaUSORIK) = TAU(K)*»2 

PRMT 201 

120 

uONT INuE 

PRMT 202 

C 


PRMT 2u3 

C 

lVALUATE FULY-'JOMI als 

PRMT 204 

C 


PRMT 2U5 


CALL PJLY(f^AX) 

PRMT 2u6 

C 


PRMT 2U7 


uQ 200 K2 = l,((.!A 

PRMT 208 

C 


PRMT 209 

c 

utTEKMINL Y(A) 

PRMT 210 

c 


PRMT 211 


CALL Y:>U6A (M(K2) ) 

PRMT 212 

c 


PRMT 2*.3 

c 

CALCULATE THtURETICAL LOG TiMtb ANO TIMES 

PRMT 214 

c 


PRMT 215 


CALL YTh(M(K2M 

PRMT 2l6 

c 


PRMT 217 

c 

COMPUTE RESIDUAL 

PRMT 2ih 

L 


PRMT 219 


CALL RlSID(M(K2) ) 

PRMT 2^:0 

c 


PRMT 221 

L 

make O.'JE entry in SUMMARY TABLE 

PRMT 22^ 

C 


PRMT 223 


ntry=;niiry+i 

PRMT 224 


TABLE! i,NTKY)=g{K5) 

PRMT 225 


TABLE! 2,NTRY)=R(K4) 

PRMT 226 


I TBLE! 3,NTKY ) =M( K2) 

PRMT Zc7 


TaBLE(4,MTRY)=TTA( K3) 

PRMT 2^8 


I ABLE! 3,NTkY)=Ya 

PRMT 229 


TABLE! 6,NTRY)=L)D 

PRMT 230 


rRG0R2=NTRY.E0.2»Ll.viES 

PRMT 231 


IF ! TKuCRP ) GJ TO 1 TO 

PRMT 23^ 


bC TU 190 

PRMT 233 

c 


PRMT 23^ 

c 

ULTPUT:> ONE PAGE OF SUMMARY TABLE 

PRMT 235 

c 


PRMT 236 

c 


PRMT 23^ 

c 

uUTPUT TITLE AND HEADINGS FOR SUMMARY TAlLE 

PRMT 238 

c 


PRMT 239 

170 

^RITE ! 6, 9005) !TITLE!K),K=1,12) , DEGREE 

PRMl 240 


IF (LINFS.EQ.51) WRITE (6,9006) K T A , KR , KG , KM , \ 

PRMT 


WRITE (6,9007) 

PRMT 242 


TRGGR1=MRY.LE.L iNEb 

PRMT 243 


LIMIT=LINES 

PRMl 244 


IF (TRuGRi) LIMIT=NTRY 

PRMT 245 


jO 180 K=l, LIMIT 

PRMT 246 


WRITE ( 6,9 008) (TABLE! I, K), 1 = 1,2), 1 TBLE! 3, K) ,!TABLE(I,K), 1-^,6) 

PRMT 247 


IF ( TRuGRl ) GU TU 180 

PRMT 246 


KUL2=K+L INES 

PRMT ^49 


IF ! TRGGR2 ) GO TU 1 75 

PRMT 230 


IF (KUL2.GT.NIRY) GO TO 180 

PRMT 23 1 

175 

WRITE (6,9009) ! TABLE! I ,K0L2) , 1 = 1, 2) , I ToLE! 3,KUL2 ) , 

PRMT 252 

1 ( TAbLt { I ,K0L2 ) , 1=4,6) 

PRMT 233 

180 

CONT 1 NUE 

PRMT 254 


ntry=o 

PRMT 255 


L INES = 55 

PRMT 256 


20 



non noo oor-oo ooroo oco 


190 


200 

300 

400 

500 


1000 

1010 

1020 

1030 


1031 

1032 

1034 

1035 


1040 


IF (TR0GR3) GU TO 1000 


prmt 

2s 7 



PRMT 

258 

SAVE VALUES PRODUCING SMALLEST 

RESIDUAL 

PRMT 

259 



PRMT 

260 

IF (DOi.LE.OU) GU TO 200 


PRMT 

2oi 

Mi = M(K2) 


PRMI 

2u2 

TTA1=TTAI K3) 


PKMT 

26 3 

R1 = R(K4) 


PRMT 

264 

U1 = 01K5) 


PRMT 

265 

yAi = YA 


PRMT 

266 

U01=DU 


PRMT 

267 

CGNT li\UE 


PKMT 

268 

CONTINUE 


PRMT 

2 o 9 

CONTINUE 


PKMl 

2 ^0 

CONTINUE 


PRMT 

2/1 

TRCGR3=.TRUE. 


PRMT 

2 

IF (NTKY.NE.O) GU TO 170 


PRMT 

? 



PRMT 

274 

END MAJOR LOUP 


PKMT 

2 Ib 



PRMT 

2 It 

OUTPUT OPTIMUM VALUES AND HEADING FOR FULL TAtiLE 

PKMT 

27 7 



PRMT 

278 

uONTiNUE 


PRMT 

2 /9 

.JRITL (6,90C:>) ( 1 1 TL t ( K. ) , K = 1 , 1 2 ) , 'JfcGKEE 

PRMT 

2F0 

L INFS=3 


PKMl 

281 

WRITE (6,90iU) CltRlf M1,TTA1, YAlfODl 


PRMT 

26^ 

L lNtS = L INES + S 


PKMT 

283 

rtRlIL (6,9011) 


PRMT 

284 

L 1NES = L lNES + 3 


PR^. I 

285 



PRMT 

286 

CALCULATE THEJRETICAL TIMES, RATIOS UF DIFFERENCE^ 

PRMT 

2b / 

TO RESIDUAL, AND VALUES uF THE 

parameter. Fur thf 

PRMT 

2 8 8 

PAKAMEIER PRODUCING THE MINIMUM 

residual 

PRMT 

289 



PRMT 

2 90 

DC 10)5 K=1,N 


PRMT 

29 1 

T01FF=a 6S( TT(K)-TTAl) 


PRMT 

292 

:>1GU(K)=SIGMA(K)**U1 


PR^ 1 

29 3 

IF (TDIFF) 1032, 1031 , 1032 


PRMT 

29A 

TAU(K) =0. 


PRMT 

2 V 5 

uG TO 1034 


PRMT 

296 

TAU(K)=SIGC(K)*TDIFF**K1 


PRMT 

2 9 / 

TAUSCK(K) = 1AU(K)**2 


PKMT 

29P 

CONT INUF 


PR^.T 

^ 19 

DD=DDi 


PRMT 

lOC 

CALL PuLY(Ml) 


PRMT 

ic 1 

CALL YsUbA(Ml) 


PRMT 

3C2 

CALL Y1F. (Ml) 


PRMT 

)0 3 

CALL RATIGI 


PRMT 

3C4 

CALL PARAM 


PRMT 

)u5 



PRMT 

) U 6 

OUTPUT FULL TABLE 


PRMT 

)C7 



PRMT 

3u8 

o 

II 


PRMT 

309 

K = K + i 


PRMT 

UO 

WRITE 16,9012) TT( K) , SIOMAIK ) , X( K) , T(K) , T IME (K) , Y( K) , YY( K) , 

PRMT 

311 

RATIU(K) ,PP(K) 


PRMT 

312 

L 1NES = L INES + i 


PRMT 

313 

IF (K.EQ.N) UU TO 1 


PRMT 

314 

IF (LINES.LT.60) GU TO 1040 


prmt 

315 

WRITE 16,9005) ( T 1 T LE ( KKK ) , KKK= 1 , 1 2 ) , 

DEGREE 

PRMT 

316 

WRITE (6,9011) 


PRMT 

31 / 

L lNES = b 


PRMT 

316 

uO TO 1040 


PRMT 

319 



PRMT 

320 

END UF PROGRAM 


PRMT 

32 1 



PRMT 

U2 


21 



ooo ooo ooo noo onooo or. non 


PKMT 

FORMAT STATFMENT3 FOR PROGRAM PRMT 

PRMl 

FORMATS FOR iiMPUT PRMT 

PRM7 

9001 FORMAT (12A6) PRMT 

9002 FORMAT (A3,F10.0) PRMT 

9003 FORMAT (A3, 12) PRMT 

9004 FORMAT ( A 3 , OPF 1 0 • 0 , 3P F 1 0 . 0 , OPF I 0 . 0 ) PRMl 

PRMT 

FORMAT^ FOR OUTPUT PRMT 

PRMT 

TITLE (SKIP3 TO NEW PAGE) PRMT 

PRMT 

9005 FOKMATl 1H1,20X, 12A6/1H ,30X,A6,10H PARAMLTEK/lH ) PRMT 

PRMT 

:iUMMARV OF INPUT PRMT 


PRMT 

9006 FORMAT (IN , 1 0 X , 4 5 HUR E t P / RU P T UR E PARAMETERS ARE I N Vt S T I l»AT ED FUR/ PRMT 
IIH ,I2,18H VaLUE(S) of T(A),,I3,25H TEMPERATURE E XP UNE N T ( S ) , , 13 , PRMT 
224H STRESS EXPONENT { S ) , ANU,I3,14H POL YN JM I A L ( S ) / 1 M , 1 0 X , 5 MU S I \jG , PRMT 


314, 12H DATA POINTS/iM ) PRMT 

PRMT 

hEAOlNbS FOR SUMMARY TABLE, UNE LINE OF SUMMARY TABLE PRMT 

PRMT 

format (1H , 2(2X, i HU, 7X,1 hr, 6X,1HM,3X,4HT{A),5X,^HY(A),^X, PRM T 

1 8HSTU,0EV. , iOX ) / IH ) PRMT 

FORMAT ( IH ,0PF5.2, F8.2, 15,F9.0,F10.2,lP£li.2) PRMT 

FORMAT (1H+,S8X,0PF5.2,F8.2,I5,F9.0,F1C.2,1PE11.2) PRMT 

PRMT 

OPTIMUM VALULo PRMT 

PRMT 

FORMAT! IH 10X44HVALUES PRODUCING SMALLEST STANDARD DE V I AT lUN/ 3H0u = PRMT 
1FS.2,4H, R = Fi>.2,4H, M=I2,7H, T ( A ) = F 6 . 0 , 7 h , Y ( A ) = F 9 . 3 , 1 1 H , S T i) . UE v . PRM 1 
2=lPt9,2/lH0) PRMT 

PRMT 

hEADINbS FOR FULL TABLE, ONE LINl OF FULL TABLE PRMT 

PRMT 

9011 FORMAT (5H T tMP , 4 X , 6H 5 TR E S S , 3 X , 3HL OG , 6 X , 4H T I M E , 6 X , 6HC AL U D , 3 X , PRMT 

13HL0G, 3X, 8HCALC L OG , 2 X6H JE V/ Su , 3 X , 9HP A R AM E T L R / i H , b X , 6 H ( * E -3 ) , 2 X , PRMT 
26HSTRES5,14X,4HTIME,5X,4 hTIME,4X,4HT1ML-/1H ) PKM1 

9012 FORMAT (IH , 0 P F 5 . 0 , E8 . 1 , F 8 . 3 , 2 F 1 C . 1 , 3 F 9 . 3 , I P F I 2 . 3 ) PRMT 

C PRMT 

9999 FORMAT ( IHI ) PRMT 

C PRMT 

END PRMT 


3^:3 
3^4 
32 3 
3^6 
3^7 
32 8 

32 9 

330 

331 

332 

333 

334 

333 
336 
33/ 

338 

339 

340 

341 

342 

343 

344 

345 
345 
347 
34fi 
349 

330 

331 

332 

33 3 

334 
353 
336 
33 / 

338 

339 

360 

361 
3 6rl 
36 3 
364 
i 6 3 
3t>6 
367 
3bE 
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$IBFI 

C 

I'C PRMtJLK L 1ST, REF, DECK 

SETS FIRST POLYNOMIAL TO UNITY AT ALL STATIONS AND STORES 

PRMB 

1 

C 

ALPHAMERIC CODE WORDS 

PRMB 

2 

C 


PRMB 

3 


BLOCK DATA 

PRMB 

A 


COMMON / PL YNML/QQ( 21,2001 , OTHERS! 85 » /ENO/END/NAMES/NAMES! 2 I 

PRMB 

5 


DATA (0Q( l.K) ,K=1,200)/200*1./,END/3HENO/, 

PRMB 

6 


1 (NAMES(K) ,K=l,2)/12HLARSONLINEAR/ 

PRMB 

7 


END 

PRMB 

8 


$IBFTC PARAM L I ST, REF , DECK 



c 

SUBROUTINE FOR EVALUATING THE PARAMETER AT EACH POINT 

PARM 

1 

c 


PARM 

2 


SUBROUTINE PARAM 

PARM 

3 

c 


PARM 

h 


COMMON /FDATA/SIGQ(200) ,TAU(200) , OTHERS (600) ,Y(200) 

PARM 

6 


1 /CALC/PP(2 00) t OTHER 1(600 )/N/N 

PARM 

6 


2 /PLYNML/0TH6R2(4221 ) f YA» 0THER3 ( 63 ) 

PARM 

7 

c 


PARM 

8 


00 10 K=1,N 

PARM 

9 


PP(K) = ( Y(K)-SIGQ(K)*YA)/TAU(K) 

PARM 

10 

10 

CONTINUE 

PARM 

11 


RETURN 

PARM 

12 


END 

PARM 

13 


ilBFTC YTH 

L 1ST, REF, DECK 



C 

SUBROUTINE FUR CALCULATING TIMES AND LOG TIMES 

FROM THE PARAMETER YTH 

1 

C 



YTH 

2 


SUBROUTINE YTH(M) 

YTH 

3 

C 



YTH 

4 


CUMMON /CALC/OTHERS (400) .TIME (200) , YY(200) 

YTH 

6 


1 

/FDATA/S IGQ( 200) , TAU(200) , OTHER 1 ( 800) 

YTH 

b 


2 

/PLYNML/QQ( 21, 200) , U ( 2 1 ) , YA , UTHER2 ( 63 ) 

YTH 

7 


3 

/N/N 

YTH 

8 

C 



YTH 

9 



uO 10 K=1,N 

YTH 

10 


YY(K) 

= 0. 

YTH 

11 

10 


CONTINUE 

YTH 

12 


Ml = 

M + 1 

YTH 

13 



UO 30 K=l,N 

YTH 

14 



UO 20 J=1,M1 

YTH 

15 


YY(K) 

= YY( K)+QQ( J ,K) *U( J ) 

YTH 

16 

20 


CONTINUE 

YTH 

17 


YY( K) 

= TAU(K)*YY(K)+SIGQ(K)*YA 

YTH 

18 


TIMF(K) = 10.»*YY(K) 

YTH 

19 

30 


CONTINUE 

YTH 

20 



RETURN 

YTH 

21 



END 

YTH 

22 


23 



non 


$IBFTC RAT 101 L I ST, REF, DECK 

SUBROUTINE FUR CALCULATING RATIOS 

OF INDIVIDUAL RESIDUALS TO ROOT-MEAN-SQUARE RESIDUAL 

SUBROUTINE RATIOl 
C 

COMMON /FDATA/OTHERS( 1000) , Y(200) 

1 /CALC/OTHERl ( 200 ), RAT 10(200) ,OTHER2(200) , YYI 200) 

2 /N/N/DD/DD 
C 

DO 10 K=1,N 

RATIO(K) = ABS( Y(K)-YY(K) )/0D 
10 CONTINUE 

RETURN 

END 


RATO 1 
RATO 2 
RATO 3 
RATO A 
RATO S 
RATO 6 
RATO 7 
RATO 8 
RATO 9 
RATO 10 
RATO 11 
RATO 12 
RATO 13 
RATO 14 


$ IBFTC RES ID L I S T , R E F , DECK 


c 

SUBROUTINE FUR CALCULATING RESIDUAL 



RESU 

1 

c 





RESU 

2 

c 

THE 

RESIDUAL IS BASED ON THE LOG OF THE 

TIME. 


RESU 

3 

c 

IT 

IS defined AS THE SQUARE ROUT OF THE 

SUM OF THE 

SQUARES 

OF RESU 

4 

c 

THE 

INDIVIDUAL RESIDUALS DIVIDED BY THE 

DIFFERENCE 

BETWEEN 

THE NUMRESU 

S 

c 

BER 

OF DATA POINTS AND THE DEGREES OE FREEDOM. THE 

DEGREES 

OF RESU 

6 

c 

FREEDOM, KK, DEPENDS ON THE PARAMETER (SEE MAIN BODY OF REPORT). RESU 

7 

L 


KK=2 FOR LARSUN-MILLER PARAMETER 



RESU 

8 

c 


KK=3 FUR LINEAR PARAMETER 



RESU 

9 

c 


KK=5 FUR GENERAL PARAMETER 



RESO 

10 

c 





RESU 

11 

c 

DD 

= SQRT( ( Y-YY)**2/(N-M-KK) ) 



RESU 

12 

c 





RESU 

13 


SUBROUTINE RESIO(M) 



RESO 

14 

c 





RESD 

IS 


COMMON /FDATA/OTHERSl 1000) , Y( 200) 



RESD 

16 


1 

/CALC/OTHERU 600 ) , YY( 200) 



RESU 

17 


2 

/DD/Du/N/N/DEGR£E/DEGREE/NAMES/FAMES( 2) 


RESU 

18 

c 





RESD 

19 



IF ( DEGREE. EQ.FAMES( 2 ) ) GO TO 20 



RESD 

2 0 



IF ( DEGREE. EQ. FAMES! 1 ) ) GO TO 10 



RESD 

21 


KK 

= 5 



RESU 

22 



GO TO 30 



RESU 

23 

10 

KK 

= 2 



RESU 

^4 



GO TO 30 



RESU 

25 

20 

KK 

= 3 



RESU 

26 

30 

D = 

N-M-KK 



RESO 

27 


DO 

= 0. 



RESU 

28 



DO 40 K=1,N 



RESD 

29 


DD 

= DD+{ Y(K)-YY{K) )**2 



RESU 

3C 

40 


CONTINUE 



RESU 

31 


DD 

= SqRT(DD/D) 



RESU 

32 



RETURN 



RESD 

33 



END 



RESD 

34 


24 


o o n o o 


SIBFTC YbUBA L I S T , R h F , 0 hCK 

SUBRUOIINE FUR EVALUATliMG Y(A) 

THIS SuBROUTI.NE ALSu EVALUATES THE UUAhTITIES» U» .MECESoaRY 
FOR DETERMINING THE THEORETICAL LOG TIMES. 

:>Ui3RUUr INE Y:>UBA(M) 

C 

COMMON /PLYNML/UQ( ^JL,2 00),U(2U,YA,A(21),a(21),C(2l) 

1 /EOATA/SIGU(200) »TAU( 200) , TAUS0R(2O0) , 0 T ME K S ( 40 0 ) , 

2 Y( 200)/N/N 



AO = 

0. 


CO = 

0. 

DO 10 K=1,N 


AO = 

AO+SlGw(K) »*2 


CO = 

Cu^-SlGUlKIi^YIK) 

10 


GONT IiNJE 


Ml = 

M + l 

DO 20 J=l,Ml 


A( J ) 

= 0. 


a( J ) 

= 0. 


C( J) 

= 0. 

20 


gQNT iNut 
uO 40 J=l,Mi 
DU 30 K=1,N 


AI J ) 

= AI J ) +J>IGU( K ) «^TAu( lO *00 ( Jf K I 


6( J ) 

- H(J)^-rAUS0R(K)*D'w(J»K)*»2 


C( J) 

= CU ) + TAU( K ) *Y( K ) * J ,K ) 

30 


gONT inue 

40 


gONT INuE 


SUMl 

= 0. 


SUM2 

= 0. 

jU 50 J=I,M1 


AOB 

= A{J)/i3(J> 


:>UM1 

= SUMH-AOa*C t J ) 


:»UM2 

= SUM2^-A0a*A M ) 

50 


1.UNT INuE 


YA = 

(C0-SUMl)/(A0-SuM2) 
DO 60 J=l,Ml 


U( J) 

= (CU)-A( JI*YA)/u( J) 

60 


CUNT INUE 

RETURN 

tND 


YSUI3 

1 

YSU6 

2 

YSUG 

3 

YSU8 

4 

YSUtt 

5 

Y:>U6 

6 

Ysua 

f 

Ysua 

6 

YSUb 


YSUb 

1C 

YSua 

11 

YSUb 

12 

YSU6 

13 

YSUB 

14 

YSU6 

16 

YSUB 

1 6 

YSUri 

17 

YSUb 

18 

YSua 

19 

YSua 

c. J 

YSUtJ 

21 

YSUG 

2.^ 

Y6U6 

^3 

YSUB 

4 

YSUB 

25 

YuUb 

t. 6 

YSUa 

27 

YSUB 

28 

YSUb 

2 9 

YSUb 

30 

YSUb 

31 

YSUB 

32 

YSUB 

33 

YSUB 

34 

YSUB 

35 

YSUB 

36 

YSUB 

3/ 

YSUB 

3B 

YSUB 

39 

YSUB 

A j 

YSUB 

41 

YSUB 

42 

Y:>UB 

4 3 
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SIBFTC PULY LIST, REF, OhCK 


c 

SUBROUTINE FUR EVALUATING ORTHObUNAL POiyNUMIALS 

POLY 

1 

c 


POLY 

2 

c 

ALL POLYNOMIALS UP TO MAXIMUM DESIRED UEGREE ARt EVALUATED 

POLY 

3 

c 

AT EACH DATA POINT 

POLY 

4 

c 


POLY 

5 

c 

THE FIRST POLYNOMIAL IS IDENTICALLY EQUAL TO UNITY 

POLY 

6 

c 

THESE VALUES ARE STORED BY A BLOCK DATA SUBROUTINE 

POLY 

7 

c 


POLY 

8 


SUBROUTINE POLYIM) 

POLY 


c 


PULY 

10 


CUMMQN /FDATA/0THERU400) ,TAUSQR(200) ,UTHER2(200) ,XX(200) , 

POLY 

1 i 


1 UTHER3(200) 

PULY 

12 


2 /PLYNML/UQ(21,200) , OTHERS! 45) , AL P HA ( 2 0 ) t uE T A { 2 0 ) 

POLY 

13 


3 /N/N 

POLY 

lA 

c 


POLY 

15 


SI = 0. 

PULY 

16 


S2 = 0. 

POLY 

17 


DO 10 K=1,N 

POLY 

16 


:»! = Si + XX(K)*TAUSOR(K) 

PULY 

1 -j 


S2 = S2+TAUSwR(K) 

POLY 

c C 

10 

CONTINUE 

PULY 

21 


ALPHAU) = S1/S2 

PULY 

22 


DO 20 K=UN 

PULY 

^3 


UQ(2,K) = XX(K)-ALPHA( L) 

POLY 

24 

20 

CONTINUE 

PULY 

25 


IF (M.lE.U RlTURN 

PuLY 

^ 6 


DO 50 K=2,M 

PULY 

dl 


SI = 0. 

POLY 

ch 


S2 = 0. 

PULY 



S3 = 0. 

POLY 

3G 


II 

o 

• 

POLY 

31 


DO 30 J=1,N 

PULY 

32 


Ul = TaUSQR! J)*QQ(K, J) 

POLY 

33 


D2 = Di*QQ(K,J) 

PULY 

34 


SI = Si+XX(J)*D2 

POLY 

35 


S2 = S2+D2 

PULY 

36 


S3 = S3+XX( J)*Dl*Qw(K-i, J) 

POLY 

37 


S4 = S4 + TAUSCR( J ) *Qu( K-1 , J ) **^2 

POLY 

38 

30 

CONTINUE 

PULY 

39 


ALPHA(K) = S1/S2 

PULY 

AC 


BETA(K) = S3/S4 

PULY 

A1 


DO 40 J=I,N 

POLY 

42 


QU(K+1, J) = (XX{ J )-ALPHA(K) )*QQ(K,J)-BETA(K) *QC(K-l , J) 

PULY 

43 

40 

CONT INUE 

PULY 

44 

50 

CONTINUE 

POLY 

45 


RETURN 

POLY 

46 


END 

POLY 

a7 
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O O O O 4A 


IBFTC SCALfc L lST,KtF,oeCK 

SUBROUTINE FUR SCALING LOGS OF STRESS 

THt SCALED VALUES LIE IN THE REGION -2 TO 2 

SUBROUTINE SuALE 

COMMON /FDATA/UTHERU 600) , X( 200) , XX( 200) »0THER2 ( 2 00) /N/N 

BIG = 0 • 

SMALL = I.E5 

00 10 K=1,N 

BIG = AMAXl(blG,X(K) ) 
oMALL = AMINK SMALL, X(K) ) 

10 CONTINUE 

A = 4./(BIG-5MALL) 

B = 2.»( b I SMALL )/( BIG-SMALL) 

00 20 K=l,N 
XX(K) = A»X{k)-B 
20 CONTINUE 

RETURN 
END 


SCAL 1 
SCAL 2 
SCAL 3 
SCAL 4 
SCAL S 
SCAL 6 
SCAL 7 
SCAL B 
SCAL 9 
SCAL 10 
SCAL ll 
SCAL U 
SCAL 13 
SCAL 14 
SCAL 15 
SCAL lb 
SCAL 17 
SCAL 18 
SCAL 19 
SCAL 20 
SCAL 21 
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1 

Index^ 

i 


TABLE I. - CALCULATION OF POLYNOMIALS FOR THEORETICAL DATA FOR THIRD DEGREE POLYNOMIAL 
[Temperature intercept, Tg^, 600*^ F. ] 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Tempera- 

Time, 

Stress, 

log t 

log cr 


Scaled 


Polynomial 


ture. 






log a. 














Si 

Op 

hr 

psi 




X 


^2 

% 


1100 

4954. 68 

56 300 

3. 69501 

4. 75051 

500 

2.0 

1 

1. 3619 

-0.19594 

-0.50845 

1100 

11365. 9 

19 800 

4. 05560 

4. 29666 

500 

1. 3806 

1 

. 30341 

-1. 6875 

-.57205 

1200 

625. 342 

30 300 

2. 79612 

4. 48144 

600 

1, 6328 

1 

.85576 

-1. 4583 

1. 7195 

1200 

2908. 

5 080 

3. 46359 

3. 70586 

600 

.57433 

1 

-.80869 

1. 5741 

-1. 2980 

1300 

117. 371 

12 900 

2.06956 

4. 11059 

700 

1.1267 

1 

. 18925 

2.5067 

-1.0745 

1300 

1340. 

778 

3. 12710 

2. 89098 

700 

-.53777 

1 

-2. 2709 

-. 90880 

. 92564 

1400 

34.4856 

4 190 

1.53764 

3. 62221 

800 

. 46017 

1 

2. 8030 

. 015445 

-. 77354 

1400 

995. 25 

66 

2. 99793 

1. 81954 

800 

-2.0 

1 

.51879 

-. 78251 

-. 98133 


TABLE II. - INTERMEDIATE CALCULATIONS FOR THEORETICAL 
DATA FOR THIRD DEGREE POLYNOMIAL 
[Temperature intercept, T^, 600® F. ] 


1 

2 

3 

4 

5 

6 

7 

Index, 

J 

a 

P 

a 

b 

c 

u 




8.0 


23. 743 


0 



3. 48 XlO® 


1 

0.27092 

0, 

5200. 

14897. 

-9. 9146X10"^ 

2 

-.57813 

1. 6548 

786. 18 

5. 7589 

2616. 

-8. 4266X10"^ 

3 

. 28432 ' 

1. 3315 

465.68 

7. 6678 

3796. 9 

-8. 1771X10"^ 

4 

. 41260 

. 80618 

286. 01 

6. 1816 

2694. 6 

-3. 6513X10-6 


29 



TABLE III. 


PIT FOR SEVERAL VALUES OF LINEAR 


PARAMETER FOR THEORETICAL DATA 


Degree of 
polynomial 

Temperature^ 

Ta 

VariaBle; 

ya 

Deviation 

2 

500 

10.54 

0. 008049 

3 

500 

10. 54 

. 009786 

4 

500 

10. 55 

. 010660 

2 

600 

9. 49 

. 004859 

3 

600 

9. 50 

. 000002 

4 

600 

9. 50 

. 000003 

2 

700 

8. 44 

. 015412 

3 

700 

8.46 

. 013937 

4 

700 

8. 45 

. 014469 


TABLE IV. - COMPOSITION OF STEELS RECEIVED 
FROM GERMAN COOPERATIVE LONG- 
TIME CREEP PROGRAM 

[As-received^ 20-mni-diani. Bar stock. ] 


Element 

Composition^ percent 



Steel 



C 

P 

K 


(23b CK) 

(I4a PA) 

(27b KK) 

Carbon 

0.065 

0. 270 

0. 068 

Silicon 

. 47 

. 26 

. 45 

Manganese 

. 60 

. 60 

. 73 

Chromium 

17.24 

2. 62 

16.14 

Molybdenum 

2.08 

. 27 

2.10 

Columbium and 

. 02 

Trace 

. 44 

tantalum 




Nickel 

11. 90 

.14 

13.12 

Titanium 

. 39 

Trace 

Trace 

Vanadium 

.10 

. 26 

.05 

Tungsten 

Less than 

Trace 

Trace 


0.005 






TABLE V. - NASA KUFTURE DATA 


(a) Steel K (27h KK) 


Temperature, 

Stress, 

Time, 

Temperature, 

Stress, 

Time, 

Temperature, 

Stress, 

Time, 



t. 

T, 

0, 

t, 


a. 


Op 

psi 

hr 

Op 

psi 

hr 

Op 

psi 

hr 

1022.00 

77 000.000 

1. 500 

1600. 00 

20 000.000 

0.400 

^'1112. 00 

60 000.000 

12. 900 

^1022.00 

72 500.000 

13. 800 

1560. 00 

20 000.000 

1.900 

®-1110. 00 

60 000,000 

34. 

^1022. 00 

72 000.000 

10. 

1520. 00 

20 000.000 

4.450 

®-1080. 00 

60 000.000 

52. 200 

^-102 2. 00 

70 000.000 

36. 700 

®’1480. 00 

20 000.000 

23. 700 

^1080. 00 

60 000.000 

37. 400 

^•1022. 00 

68 000,000 

60. 400 

^460. 00 

20 000.000 

25.500 

^1050. 00 

60 000.000 

239. 

^1022.00 

66 000.000 

73. 300 

®1440. 00 

20 000.000 

38. 

^1030. 00 

60 000.000 

445. 

^1022.00 

66 000.000 

107. 600 

®1400. 00 

20 000.000 

136. 000 

^1022. 00 

60 000.000 

989. 900 

®-1022. 00 

65 000.000 

201. 300 

a-1360. 00 

20 000.000 

394. 800 

^1020. 00 

60 000.000 

817.500 

^1022.00 

62 500.000 

250. 400 

^1340.00 

20 000.000 

704. 600 

1040. 00 

75 000.000 

.330 

®-1022. 00 

60 000.000 

990. 

^1320. 00 

20 000.000 

1212. 

1022.00 

75 000.000 

5.850 

®1022, 00 

60 000.000 

817.500 

1320. 00 

40 000.000 

2. 700 

^■1000. 00 

75 000.000 

15. 600 

^1022.00 

55 000.000 

3 680. 

1290. 00 

40 000.000 

7.500 

^980. 00 

75 000.000 

46.500 

1112. 00 

68 000.000 

.750 

^260.00 

40 000.000 

15. 200 

^■960. 00 

75 000.000 

138. 

1112, 00 

65 000.000 

2.250 

^■1230. 00 

40 000.000 

44. 400 

^940. 00 

75 000.000 

542. 

1112. 00 

62 500,000 

4. 300 

^1170. 00 

40 000.000 

377, 

®-920. 00 

75 000,000 

579. 600 

^1112. 00 

60 000.000 

13. 900 

^1140. 00 

40 000.000 

1417. 

^1120. 00 

50 000.000 

186.100 

^1112. 00 

57 500.000 

22. 700 

^1125.00 

40 000.000 

2110. 

^1200. 00 

40 000.000 

130. 200 

^1112. 00 

55 000.000 

51, 500 

^ 112. 00 

40 000.000 

5367. 

^1280. 00 

30 000.000 

132. 700 

®-1112.00 

52 500.000 

147,500 

1200. 00 

60 000,000 

.610 

^1340, 00 

25 000.000 

125.800 

112. 00 

50 000.000 

2R3. 

1170. 00 

60 000.000 

1.250 

^1500.00 

15 000.000 

51. 300 

^1112. 00 

45 000.000 

1 020. 

1150. 00 

60 000.000 

4. 400 

^1560. 00 

12 000.000 

41, 700 

^1112.00 

43 000.000 

1 579. 

1140. 00 

60 000.000 

4.500 

^1580.00 

10 000,000 

32. 400 

1112. 00 

37 000.000 

15 140, 

^1120.00 

60 000.000 

10. ?00 

^1540.00 

10 000.000 

148. 200 


(b) Steel G (23b CK) 


Temperature, 

Stress , 

Time, 

Temperature , 

Stress, 

Time, 

Temperature, 

Stress, 

Time, 

T. 

cr. 

t j 

T, 

a. 


T, 


t, 

Op 

psi 

hr 

Op 

psi 

hr 

Op 

psi 

hr 

^1600. 00 

5 000.000 

570. 200 

^1230.00 

30 000.000 

175. 700 

®-1202.00 

36 000.000 

68. GOO 

^1620. 00 

5 000.000 

186. 600 

^1250. 00 

30 000.000 

103, 500 

^1202.00 

38 000.000 

59. 300 

^1660. 00 

5 000.000 

156. 800 

ai280.00 

30 000.000 

58. 100 

^1202.00 

42 000.000 

24. 400 

^1680. 00 

5 000.000 

91. 600 

1292.00 

30 000.000 

21 300. 

^'1202. 00 

44 000.000 

14.500 

^1700. 00 

5 000.000 

62. 700 

^1310. 00 

30 000.000 

22. 500 

^1202.00 

45 000.000 

22. 900 

^1740. 00 

5 000.000 

40. 500 

^1112.00 

40 000.000 

667. 900 

1202.00 

46 000.000 

7. 

^•1780.00 

5 000.000 

10. 600 

^•1120. 00 

40 000.000 

785. 400 

1202.00 

48 000.000 

2.850 

^1425.00 

10 000.000 

1690. 

^1150.00 

40 000.000 

266. 700 

1202.00 

49 000.000 

2.550 

^1450. 00 

10 000.000 

550. 300 

^1170. 00 

40 000.000 

127. 800 

1202.00 

50 000.000 

1.470 

^•1480. 00 

10 000.000 

270. 

^1202. 00 

40 000.000 

44. 100 

8 

00 

1 — 1 
a 

18 000.000 

859. 700 

^1500. 00 

10 000.000 

170. 

^1202. 00 

40 000.000 

74. 

^1292.00 

23 000.000 

194. 600 

^1520.00 

10 000.000 

128.500 

^1210.00 

40 000.000 

40. 500 

^1292.00 

25 000.000 

75. 

^1560. 00 

10 000.000 

40. 

^1220. 00 

40 000.000 

37. 800 

^12'.^2.00 

28 000.000 

34 . GOO 

^1570.00 

10 000.000 

31. 500 

^1 240. 00 

40 000.000 

17. 200 

^1292.00 

29 000.000 

31. 

^1600. 00 

10 000.000 

15.800 

1270.00 

40 000-000 

4.500 

^1292.00 

32 000.000 

13. 300 

1650. 00 

10 000.000 

5. 250 

1280.00 

40 000.000 

1. 200 

■^1292.00 

33 000.000 

19. ROO 

1700. 00 

10 000.000 

1. 750 

1292.00 

40 000.000 

1. 300 

^1292.00 

34 000.000 

10. 400 

^1202. 00 

20 000.000 

3307. 

1300. 00 

40 000.000 

.800 

1292.00 

36 000.000 

2.750 

^1260. 00 

20 000.000 

667. 400 

^1112.00 

34 000.000 

2 274. 

1292.00 

37 000.000 

7.600 

^1290. 00 

20 000.000 

255. 

^1112.00 

43 000.000 

363. 100 

1292.00 

38 000.000 

1.650 

^1 292. 00 

20 000.000 

347. 100 

^1 112. 00 

46 000.000 

233. 900 

^1060. 00 

60 000.000 

42.500 

^1292. 00 

20 000.000 

363. 

^1112. 00 

46 000.000 

261. 400 

^1300.00 

25 000.000 

89.600 

^•13 20. 00 

20 000. 000 

180, 400 

^1112. 00 

48 000.000 

183. 100 

®-1360, 00 

19 000.000 

95. 

®'1360. 00 

20 000.000 

82. 

^1112. 00 

50 000.000 

84. 500 

^■1430. 00 

15 000.000 

71. 400 

®'1400. 00 

20 000.000 

28. 900 

^1112. 00 

52 000.000 

65. 600 

®-l480. 00 

12 000.000 

147. 900 

1440. 00 

20 000.000 

9. 

^1 112. 00 

54 000.000 

39. 300 

^1570. 00 

8 000.000 

104. 

1480. 00 

20 000.000 

2.500 

^1112.00 

57 000.000 

23. 300 

^1650. 00 

6 000.000 

140. 900 

^1112.00 

30 000.000 

4258. 

^1202. 00 

25 000.000 

1 074. 

^1140. 00 

34 000.000 

1077. 

^1160. 00 

30 000.000 

1110. 

^1202. 00 

34 000.000 

109. 400 

^1320. 00 

15 000.000 

1505. 

^1180. 00 

30 000.000 

696. 300 

^•1202. 00 

35 000.000 

124. 300 

®'1480. 00 

8 000.000 

2237. 

^■1 202. 00 

30 000.000 

350. 




^1540.00 

6 000.000 

1258. 


^Data point used in parametric analysis. 



Q2ABLE V. - Concluded. NASA RUPTUEE DATA 


(c) Steel P (I4a PA) 


Temperature , 

Stress, 

Time, 

Temperature, 

Stress, 

Time, 

Temperature, 

Stress, 

T, 

Oy 


T, 

Oy 

-^y 

T, 

Oy 

Op 

psi 

hr 

Op, 

psi 

hr 

op 

psi 

932.00 

65 000.000 

3.800 

ai250. 00 

10 000.000 

19. 200 

^740. 00 

90 000.000 

^932. 00 

60.000.000 

14. 150 

Q’1220. 00 

10 000.000 

42. 

^785. 00 

80 000.000 

^932.00 

60 000.000 

14. 400 

^1180. 00 

10 000.000 

167. 

^820. 00 

70 000.000 

^932.00 

57 500.000 

10. 

^1170. 00 

10 000.000 

203. 400 

^880. 00 

60 000.000 

^932. 00 

55 000.000 

18. 900 

^'1140. 00 

10 000.000 

608. 

^932.00 

50 000.000 

^■932.00 

52 500.000 

51. 

^1090. 00 

10 000.000 

2639. 

^1022.00 

30 000.000 

^'932. 00 

40 000.000 

623. 

1100. 00 

40 000.000 

1.300 

^1050. 00 

25 000.000 

®'932.00 

30 000.000 

7 592. 

1080. 00 

40 000.000 

2. 200 

^-1090. 00 

20 000.000 

932. 00 

27 000.000 

11 410. 

1060. 00 

40 000.000 

4. 300 

^1090.00 

20 000.000 

1022.00 

58 000.000 

.580 

1050. 00 

40 000.000 

6. 800 

^1120. 00 

16 000.000 

1022. 00 

55 000.000 

. 717 

1040. 00 

40 000.000 

7.400 

^■1160. 00 

13 000.000 

1022. 00 

50 000.000 

1. 280 

^1020. 00 

40 000.000 

22. 500 

®'1230. 00 

8 000.000 

1022. 00 

47 000.000 

2.450 

^1010. 00 

40 000.000 

20. 100 

3-1290. 00 

5 000.000 

1022.00 

47 000.000 

6.200 

s-lOOO. 00 

40 000.000 

63. 300 

®'740. 00 

80 000.000 

1022. 00 

45 000.000 

3.500 

^990. 00 

40 000.000 

51. 200 

^780. 00 

70 000.000 

1022.00 

42 500.000 

6.300 

^980. 00 

40 000.000 

80. 600 

3830. 00 

60 000.000 

^1022. 00 

40 000.000 

22.500 

^960. 00 

40 000.000 

192. 100 

3880. 00 

50. 000. 000 

^1022. 00 

37 500.000 

12. 

^940. 00 

40 000.000 

427. 900 

3932.00 

35 000.000 

^1022.00 

25 000.000 

382. 200 

^930. 00 

40 000. 000 

623. 

^980. 00 

30 000.000 

1415. 00 

10 000.000 

.170 

^00. 00 

40 000.000 

2572. 

^■1000. 00 

25 000.000 

1340. 00 

10 000.000 

1. 500 

932.00 

70 000.000 

1. 400 

31030.00 

20 000.000 

1315.00 

10 000.000 

3. 700 

897.00 

70 000.000 

5. 800 

31070. 00 

16 000.000 

1290.00 

10 000.000 

6.100 

^860. 00 

70 000.000 

31. 200 

31150. 00 

8 000.000 







31220. 00 

5 000.000 


^ata point used in parametric analysis. 


Time, 

lir 

57.100 

84. 

195. 800 

120 . 

103. 500 

186. 700 

123.500 

79.500 

112. 400 

183. 500 

100. 300 
97. 900 
139. 700 
996. 600 
1122. 

948. 800 
599. 

1902. 

754. 800 
970. 700 

1084. 

804. 800 

948. 500 
960. 
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TABLE VI. 


GERMAN RUPTURE DATA 


Temperature^ 

Stress^ 

Time^ 

Temperature , 

Stress^ 

Time^ 

Si 


■t. 




Op 

psi 

hr 


psi 

hr 

steel K (27b KK) 

steel 0 (23b OK) 

1022. 00 

76 899.999 

0. 100 

1292.00 

17 800.000 

1 100. 

1022.00 

66 899.999 

160. 

1292.00 

21 400.000 

300. 

1022.00 

55 500.000 

2 000. 

1292.00 

21 400.000 

250. 

1112.00 

72 500.000 

.100 

1292.00 

27 000.000 

180. 

1112. 00 

64 000.000 

10. 

1292.00 

27 000. 000 

140. 

1112. 00 

55 500.000 

35. 

1292.00 

47 000.000 

.100 

1112. 00 

44 100.000 

2 100. 


— 


1112. 00 

28 4=00.000 

52 000. 

£ 

teel P (l4a I 

>A) 

steel C (23b C 

K) 

932.00 

84 000.000 

0.100 




932.00 

75 500. 000 

. 100 

1112.00 

14 200.000 

60 000. 

932.00 

78 399. 999 

2. 

1112. 00 

17 800.000 

30 000. 

932. 00 

55 500.000 

150. 

1112. 00 

28 400.000 

3 500. 

932.00 

44 100.000 

1 700. 

1112. 00 

28 400.000 

3 000. 




1112.00 

28 400.000 

2 200. 

932.00 

34 200.000 

2 600. 




932.00 

27 000.000 

16 000. 

1112.00 

35 600.000 

1 200. 

932. 00 

22 800.000 

22 000. 

1112. 00 

44 100.000 

520. 

932.00 

17 100.000 

100 000. 

1112.00 

51 200.000 

150. 

1022.00 

72 599.999 

.100 

1112. 00 

59 800.000 

.100 




1202. 00 

11 400.000 

82 790. 

1022.00 

69 699. 999 

.100 




1022.00 

65 500.000 

1. 200 

1202.00 

14 200.000 

15 000. 

1022.00 

59 800.000 

1. 500 

1202.00 

17 800.000 

6 500. 

1022.00 

35 600.000 

150. 

1202. 00 

22 800.000 

1 800. 

1022.00 

27 000.000 

300. 

1202. 00 

28 400.000 

550. 




1202.00 

35 600.000 

124. 

1022. 00 

22 800.000 

400. 




1022.00 

22 800.000 

900. 

1202.00 

42 700.000 

5. 

1022. 00 

17 100.000 

2 100. 

1202.00 

52 600.000 

.100 

1022. 00 

13 900.000 

6 500. 

1292.00 

, 11 400.000 

30 000. 

1022.00 

13 900.000 

8 000. 

1292.00 

11 400.000 

20 000. 

1022.00 

11 100.000 

10 000. 

1292.00 

13 900.000 

4 500. 

1022. 00 

8 830. 000 1 

68 000. 



TABLE VTI. - CREEP DATA FOR COLUMBIUM ALLOY FS-85 


Temperature ^ 

Stress^ 


Time; t; hr 


T, 





Op 

psi 

l“Percent creep 

2-Percent creep 

5 -Percent creep 

2005 

25 000 

0. 6 

3.0 

6.1 

1900 

25 000 

26. 

33. 

45. 

1790 

25 000 

210. 

257. 

332. 

2175 

18 000 

4. 9 

7. 8 

13. 

2400 

10 000 

3.4 

5. 7 

10. 8 

2300 

10 000 

25.4 

41. 

68. 

2200 

10 000 

54. 

84. 

133. 

2100 

10 000 

355. 

500. 

765. 

2100 

10 000 

380 

570. 

875. 

2000 

10 000 

775. 

1325. 

2175. 

2000 

10 000 

900. 

1420. 


2000 

8 500 

2480. 



2575 

6 000 

5.6 

10. 

22. 2 

2200 

6 000 

425. 

710. 

1370. 

2800 

4 000 

3.4 

6. 4 

13.5 

2620 

4 000 

14. 4 

26. 

56. 

2200 

4 000 

1140. 



2900 

3 000 

2. 6 

5. 4 

13. 8 

3000 

2 000 

4. 6 

9.5 

33. 2 

2450 

2 000 



950. 
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*’T!he aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena m the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concernmg its activities and the results thereof.” 

• — ^National Aeronautics and Space Act of 1958 
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